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Abstract. - We argue that the chaotic temperature effect predicted in Ising spin glasses 
should be stronger when one considers continuous (XY, Heisenberg) kind of spins, due to bigger 
entropic fluctuations. We then discuss the behavior of 3d spin glasses using the Migdal-Kadanoff 
Renormalization Group for Ising and XY spins, where we show explicitly that the chaotic length 
scale, the length beyond which equilibrium configurations are completly reshuffled when one 
changes temperature, could be far smaller for XY than for Ising spins. These results could thus 
explain why experiments are seeing a stronger rejuvenation effect for continuous spins. 



One of the most striking features of spin glasses [1] are the phenomena of rejuvenation and 
memory [2] when temperature hysteresis cycles are performed. When a spin glass approaches 
equilibrium, it ages and reduces its susceptibility; however, if one lowers the temperature after 
one such aging at temperature T±, one sees a restart (a rejuvenation) of the susceptibility 
aging at temperature T2 < T±, while memory of the previous aging can be retrieved when 
heating back. More impressive is maybe the fact that many such memory cycles can be 
performed [3]. One of the possible explanations is to consider the temperature chaos effect, 
that is the possibility that equilibrium spin configurations at temperatures T\ and T2 are 
uncorrelated, a property that has been predicted to hold in spin glasses many years ago [4, 
5]. The presence of memory is then compatible with rejuvenation because of length scale 
separation: dynamics at T2 < T\ is so slow that it is not fast enough to destroy the previous 
ordering (see for instance the ghost domains of [6]). In fact most of spin glass dynamics 
(see [7] for a classic review) could be explained if one assumes that there is a growing coherence 
length scale during aging, a strong sensitivity of equilibrium configurations to temperature and 
a separation of length scales at each temperature [8] . Another very important experimental 
result is that rejuvenation is stronger for Heisenberg than for Ising spins [3,9-11] : the sharpness 
of this effect appears to decrease continuously with the spin anisotropy. 

The purpose of this paper is to argue that this might be simply due to a stronger tempera- 
ture chaos for continuous (XY, Heisenberg) spins, and to provide a simple model to illustrate 
this assumption. Also, since it has been explained recently "why temperature chaos is hard 
to observe" [12], we would like to show how it may actually be easier to see. The paper is 
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organized as follows: first, we discuss generally temperature chaos in spin glasses, then we 
introduce a 3c? spin glass model interpolating between XY and Ising spins and the Midgal- 
Kadanoff Renormalization Group we will use to study it. We will present results showing that 
chaos is stronger when considering XY type of spins, and more generally when there are more 
degrees of freedom in the system. We will conclude by a general discussion. 

Temperature Chaos in Spin Glasses. - Let us start by the simple thermodynamic argu- 
ment of [5], and consider a 3d Edwards- Anderson [13] spin glass in the scaling/droplet [4, 14] 
approach: take two equilibrium states at temperature T < T c differing by a very large droplet 
of characteristic size I. Then the two states have free-energies that differ by AF(T) = AE — 
TAS w T (T) £ e , where T (T) is the energy stiffness coefficient. When one changes the tem- 



perature by ST then AF(T+ST) « AE- (T + ST)AS so that AF(T+ST) « T (T) l e — 5TAS. 



In this phenomenological approach, the entropy difference is associated with the droplet 's sur- 
face so that AS has a random sign and a typical magnitude a (T) £ ds / 2 7 where a (T) is called 
the entropy stiffness and d s is the fractal dimension of the droplct's surface. If d s /2 > 9, 
which follows from droplet theory, then AF(T + ST) can change sign between T and T + ST 
for length scales greater than 



with £ = d s /2 — 9. So, when temperature is changed, equilibrium configurations are changed 
on scales greater than £ c . The crucial point in the argument lies in these large cancellations 
between E and 5* in the equilibrium value of F: these make the equilibrium state very sensitive 
to changes in T. This picture of chaotic temperature dependence has been largely confirmed 
within various scaling approaches [15-17], in the naive mean field approximation of the spin 
glass phase [18] as well as in a number of toy models [18-21]. There are now analytical results 
for such a dependence in mean field models [22] but, however, weak enough to be very hard 
to see in simulations [23]. In the case of finite dimensional Ising spin glasses, it is still not 
clear if there is (or not) temperature chaos [23-25], probably because l c is quite big in that 
case, as claimed by [12]. Thus, temperature chaos does probably exist in finite dimensional 
spin glasses but is very hard to see and to confirm so that experimental relevance could be 
(and has indeed been [25,26]) questioned if one has to consider length scales bigger than those 
that could be reached experimentally. Indeed, other explanations for rejuvenation have been 
provided [8,26,27] and debates and controversies are still open [10,28]. The experimental fact 
we are interested in here is, again, that in experiments rejuvenation is found to be stronger 
for Hciscnberg than for Ising spins. 

A simple explanation would thus be that chaos is also stronger for continuous spins. Why 
should it be so? Looking at eq.(l), one sees that chaos will increase if energy (exponent 
or stiffness) decreases and if entropy (exponent or stiffness) increases. This has a simple 
interpretation: chaos is due to fluctuations of energy and entropy (remember that, in droplet 
theory, there are more fluctuations in free energies when 9 decreases); therefore the more 
fluctuations there are, the more chaos is important! One then may expect that Heisenberg 
or XY spins, which have many more degrees of freedom, will indeed have more entropy, thus 
lead to more fluctuations. The purpose of the following study is to illustrate this point in 
(some) special cases. What we need for studying temperature chaos is a spin glass model 
which is both (i) tractable enough in order to be able to compute £ c and to study the chaotic 
temperature dependence and (ii) which allows us to study continuous as well as Ising spins. 
Following early studies on spin glasses [15,16,29], a natural answer is to use the Migdal- 
Kadanoff Renormalization Group (MKRG) [30] . 




(1) 
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Fig. 1 - Left: Construction of a hierarchical lattice, for dimension d one needs 2 d_1 branches. Right: 
entropy and energy of an Ising spin glass with ±1, Gaussian and power-law tail distributions, rescaled 
to obtain T c = 1; energy decreases and entropy increases while distributions get broader, but expo- 
nents (here plotted with 6 = 0.26 and d s /2 — 1) remain the same. The inset shows the distance 
between two replicas at T = 0.1 and T — 0.11 versus L: broader distributions lead to more chaos. 

Migdal-Kadanoff Renormalization Group. - It is sometimes convenient to see the MKRG 
as an exact resolution on a hierarchical lattice (see fig.l). Since it is not possible to write 
directly an exact recursion relation for continuous spins, we will use a q-states clock model, as 
introduced by [31] for which the two limits q = 2 (Ising spins) and q = oo (XY spins) display 
a similar droplet-like ordering with a finite T c as well as the expected chaotic temperature 
dependence [17]. We will try here to be more quantitative, following the recent work [12] 
in the Ising case. Let us define the g-states clock model, following [31], with a very general 
random gauge Hamiltonian 

N 

# = -E J ^-^-) ( 2 ) 

<ij> 

where J^j(S<f>) could be any function defined for 5<f> having values restricted to 2irk/q, with 
k = 0, 1, 2, (g — 1). It is then possible to write a one dimensional decimation procedure, for 
three consecutive spins a,b,c on a chain [31] 




9-1 

F(k,T) =^exp((j a6 (0 + J bc (mod ((q + k -I), q)))/Ty (3) 

1=0 

To do a complete MKRG step, one has to pick up 2 d_1 links decimated according to (3), to 
create 2 d ~ 2 new links and merge them into the final new J, this final step being just a link 
addition. We thus have to generate a set of J^(k) at each renormalization step n; since we 
are dealing with a random system, we will do that numerically using the pool method of [15], 
the only tricky points being to take care of the numerical representation of infinitesimal and 
exponential numbers, and to check that the pool size is not too small. One can use any specific 
form for the initials conditions J^(4>); we will choose here the random gauge XY spin glass, 
thus starting with = cos(0i — 4>j — Aij), where Aij are randomly chosen between and 
2tt. Of course the initial cosine distribution is not conserved after iterations and converges to 
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a random distribution; this is indeed happening for almost any initial distribution, leading to 
universality. An exception to this rule arises when the interaction has a reflexion symetry, like 
for instance JijSiSj] in that case there is no finite T transition for d = 3, but any weak random 
perturbation changes the RG flow to that of the gauge model (see [17,31] for discussions on 
these points) which seems, therefore to be the more natural model to study. 

Energy, entropy, distance and chaos. - For a complete description of method and for- 
malism that will be used in the following, especially concerning the pool method, we refer the 
reader to [15]. The idea is to generate a representative set of J™ at each iteration step. It is 
then possible to compute all quantities of interest. The droplet free energy and entropy reads 

With q-statcs, one can choose to compute these quantities for the maximum J{4>) or to do 
an angle average; results are not significantly modified. The very standard way of studying 
temperature chaos is then to consider the length scale at which two replicas of the same system 
at temperature T and T + ST decorrelate; we will thus define l c as the length 2"° at which 
d(no) = 1/e, where the distance d(n) is defined by: 

,2, x ({J n {T)-.r{T + 5T)f) 

a (n) = 5 7T- (5) 

((J"(T)) 2 ) + ((J"(T + AT)) 2 } 

When ST « T this length scale is known to be well estimated by cq.(l) (see for example [12, 
15]), however we do not really care about the preciseness of this particular equation and 
therefore formula (5) will be used in the following to define £ c . We will now show results 
obtained from the renormalization process, first for the Ising then for the XY model. 

Ising spin glasses and the non universality of l c . In a recent and very interesting 
paper [12] l c was computed in 3d Ising spin glasses, first using the MKRG approach, then 
directly on a 3d lattice by considering interfacial energy and entropy and using eq.(l). The 
authors conclude that £ c is so big that it will be hard or almost impossible to sec it in 
numerical simulations. Let us re-examine some of the aspects of chaos in the Ising case. The 
9 exponent is well known for Ising MKRG, and even analytical studies are possible [32] giving, 
for d = 3, 9 « 0.26. The surface fractal dimension d s is also know to be d — 1 on such MK 
lattices. We have considered variations of the chaos intensity when one uses ±J, Gaussian 
and also a power-law coupling distribution (decreasing like x 1+0 5 ), and rescaled all energy 
distributions in order to obtain T c = 1. We worked at T = 0.1 and computed droplet energies 
and entropies, showing as expected that 9, as well as the entropy exponent, are the same for 
all models. The important point (see fig. 1 ) is that, as tails increase in the initial distribution, 
the energy stiffness decreases and the entropy increases. We thus expect that chaos increases 
with tails, and, computing distances, we see that this is indeed the case. Thus, in complete 
agreement with eq.(l), we find that chaos increases when the energy decreases and the entropy 
increases. More tails mean that fluctuations are stronger so, in other words, chaos increases 
with fluctuations in the system: the more fluctuations there are, the more chaos there is. 

To conclude this preliminary work, it is worth making the following comments: (a) although 
it may seem a bit trivial, we see here that £ c is not a universal quantity and indeed can change 
from more than a factor 10 even in simple 3d Ising model; it is thus somehow a bit strange 
to compare experiments and simulations in such a precise way as to consider exact values of 
length scales to rule out or confirm scenarios, (b) Another important practical information is 
that, when one wants to find chaos in a system, or when doing dynamical simulations to see 
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Fig. 2 - Left: domain wall free energy of the q-states model at T = 0.01 for different q and the 0.26 
slope. In inset T c (q) converges fastly to » 0.47. Right: interfacial entropy at T = 0.01, notice the 
non monotonic effect; convergence is reached for q > 40 . In the inset the entropy exponent converges 
to ?s 1.26 (entropy was computed using 5T — 10 -11 ). 



rejuvenation, Gaussian couplings should be used instead of discrete ± J; in dynamical studies 
especially, where the coherence length growth is very slow, the difference predicted here should 
be very important. (c) Looking at the inset of fig.l, one should also notice that the chaotic 
length is not really well defined in the sense that there is no sharp transition between £ < £ c 
and £ > £ c . Decorrelation is in fact continuous and there are also some changes for £ < t c . 
This is a crucial point (see for instance [8, 10]) since it explains why rejuvenation can appear 
before the equilibrium length scale reaches £ c , as shown explicitly in [20,21]. 

XY spin glasses and the importance of chaos. - We turn now to the q-states model and 
its XY limit. We will use the random gauge spin glass which, as opposed to the usual XY spin 
glass, has both its Ising (q = 2) and XY (q — oo) limits that display similar ordering with a 
finite T c on such lattices [17,31]. We thus start from a set of initial random link variables 

Jij(k) = cos(2irk/q — Ay). (6) 

A very nice study of chaotic properties of this model can be found in [17]. We will try here 
to be more precise, focusing on the values of £ c at the q — > oo limit. Let us first see how do 
interfacial entropy and free energy change with q. A first very important property is that the 
droplet exponent 6 fs 0.26 is the same for all q (see fig. 2) suggesting that T c is always finite; 
we have indeed checked that T c remains finite at large q. We can thus compare, changing 
from Ising to XY, similar kinds of system with a similar low temperature droplet spin glass 
phase; this will allow us to isolate the effect of continuous spins on the chaotic length. 

If the free energy has the same 9 exponent, the interfacial entropy is found to be non 
universal. However, the way the limit q — > oo is reached is a bit tricky due to a non monotonic 
effect: entropy first increases, then decreases and saturates (see fig. 2). This effect is less 
strange after the following remark: if one starts with a random distribution of </(</>) using a 
uniform (instead of cosine) distribution then the effect disappears; so this non monotonic effect 
comes from a pre-asymptotic regime, where the system is loosing its initial cosine coupling 
distribution. We need therefore to go to very large q (here q > 40) to see the XY regime. In 
contrast the entropy exponent grows monotonically from d s /2 = 1 to w 1.26, in agreement 
with the result £ « 1 for the chaos exponent when q — -> oo [17]. Finally, as expected, the 
interfacial entropy is bigger for XY spins and we can expect that so is chaos. 
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Fig. 3 - Left: distance (see eq. (5)) between equilibrium configurations at temperature T\ — 0.1 and 
T2 = 0.2 for different q. Symbols for q — 16, 32, 64 are hidden by those for q = 128, showing the 
convergence to the high q limit. In the inset the same with Ti = 0.05 and T2 = 0.06. Right: chaotic 
length £ c for the XY model and Ising models. Note that £ c can be about 10 times smaller for XY 
than for Ising spins. This plot was made using ST = 0.01 (ST — 0.1 in the inset). 



Now how does all that indeed translate in terms of the chaotic length? In [17] only the 
chaos exponent for ST — > was considered and to go beyond that we need to compute l c . In 
fig. 3 we plot the distance for different q, showing that we have again the non-monotonic effect: 
the distance first decreases, re-increases, and finally saturates with q so that great care has 
to be taken for the limit q — » 00. That said, our final results clearly show that l c is smaller 
for XY than for Ising spins (the effect would have been even stronger starting from a uniform 
distribution, thus supressing the pre-asymptotic regime) and that it can decrease by more 
than one decade. This is illustrated in fig. 3 where we plot, for many temperatures in the Ising 
and the XY limit case, l c for ST = 0.1 and ST — 0.01 (the XY limit was obtained with q = 64 
for ST = 0.1 and q = 128 for ST = 0.01). The conclusion is thus that, within this approach, 
the XY is more chaotic than the Ising spin glass, as expected from our previous argument. 

Discussion. Our main message is that, in a droplet/scaling approach, temperature 
chaos is far stronger for continuous than for Ising spins. Since such a chaotic temperature 
dependence does give rise to both rejuvenation and memory, as shown in ad-hoc models [6,33] 
as well as in dynamical MKRG approaches [20,21], one can thus use it to interpret exper- 
iments [33]. That chaos is stronger for continuous spins then may explain naturally why 
rejuvenation effects are also experimentally stronger for continuous spins, a result sometimes 
presented as a riddle. Of course, since any simple picture with a growing length scale and 
changes in equilibrium configurations with temperature are sufficient to account for that kind 
of experiments [8], there are other possible explanations for these effects [26,27], but it is fair 
to say that they do not tell anything about differences between Ising and continuous spins; this 
is therefore an indication that chaos does play a role in these memory/rejuvanation protocols. 

Note however, the very nature of the spin glass phase still being unknown and object of 
many debates, that all these scaling considerations might turn out too naive for real systems. 
We considered in this letter only the simplest possible model, but we believe nevertheless 
that, if temperature chaos does exist, its generic behavior should be well described by this 
approach. Another potential problem concerns the nature of the spin glass phase with XY or 
Heisenberg spins, which could be different from the Ising one (like the chiral ordering of [34], 
however jeopardized in most recent numerical work [35]). 
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Finally, our results also suggest some natural extensions: (i) Monte Carlo studies should 
be able to see chaos more easily in XY and Heisenberg than in Ising spin glasses, (ii) one 
may also play with the couplings to have a shorter £ c , since most work is done with ± J and 
that our results suggest instead to use a Gaussian or even a more tailed distribution, and (iii) 
dynamical studies of MK Ising spin glasses [20,21] may be usefully extended to XY spins. 
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